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Abstract. 

We have fitted the HST images and the spectral energy distributions 
(SEDs) of three proto-planetary nebulae (PPN) with a 2-D radiation 
transfer model. The geometric and mass loss properties of these PPN are 
also derived. 



1. Introduction 

The detection of bipolar reflection nebulosities in PPN suggests that the dust 
envelopes in PPN are highly asymmetric, (see Hrivnak, these proceedings). 
In order to model these objects, a 2-D radiation transfer model is necessary. 
However, a fully self-consistent determination of the source function in 2-D is 
impractical not only because of the large computing time required, but also in 
the lack of knowledge of the physical details of the scattering process. Since 
the morophology of the nebulae will be primarily determined by the geometry 
and orientation of the envelope, we have developed an approximate solution to 
simultaneously fit the SED and images of a centrally-heated dust envelope. In 
this paper, we report the models for 3 PPN, IRAS 17105-3224, IRAS 18095-^2704 
and IRAS 17441-2411, where HST images are available. 



2. The Model 



The dust envelope is assumed to be axial symmetric where the density dis- 
tribution p(r, 9) is assumed to have radial cutoffs at Tin and rout- The plane 
perpendicular to z (9 = 90°) is referred to as the "equatorial plane" and the 
directions along the z axis {9 = 0° and 180°) are referred to as the two "poles". 
The density is assumed to decrease from the equator to the poles in the form of 
a power law {p oc 9^r~'^). In order to produce the searchlight beams observed 
in the PPN, IRAS 17150-3224 (Kwok et al. 1998), cavities can be put in the 
density distribution simply assuming a open-cone structure, where p drops by a 
factor of Tscaie insidc the cone. A disk can also be put in the density distribu- 
tion in order to reprduce the dark lane observed in IRAS 17441-2411 (Su et al. 
1998). The viewing angle i is defined as 0° if the object is viewed along the pole 
(pole-on), and 90° if it is viewed along the equator (edge-on). 
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SED Plot of IRAS 17150-3324 
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SED Plot Of IRAS 18095 + 2704 
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SED Plot Of IRAS 17441-2411 
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Figure 1. The SED and model spectra for (a) IRAS 17150-3224 (b) 
IRAS 18095-^2704 (c) IRAS 17441-2411. The square and * are ground- 
based and HST photometric measurements respectively. The IRAS 
LRS is plotted as a dashed line and the model spectra as dotted lines. 



2-D model of PPN 
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We first solve the dust temperature distribution at the "polar" and "equa- 
torial" directions from 1-D radiation transfer models. The values at other angles 
are then obtained by interpolation using a power law: 

TV 

for 9 < 7r/2 
T{tt -6)= T{e) 

The three power indices 7, A^) are adjusted until best fits are obtained for 
both the SED and the images. Specifically, 7 can be constrained by the infrared 
color ([60]-[100] jim) because the total amount of dust is fixed. (3 determines the 
degree of asymmetry and is constrained by the observed width of the reflection 
nebula. is basically the same as (5, but can be adjusted in order to get the 
total energy conserved. The viewing angle i can be determined by comparing 
the flux ratio of the two reflection lobes in the simulated image to the one in the 
observed image. 

IRAS 17150-3224 and IRAS 18095+2704 are both 0-rich objects, as evi- 
denced by the presence of the 9.7/im silicate feature in the IRAS Low Resolution 
Spectra (LRS); therefore, we use silicate dust grains in the fitting. The strength 
of the silicate feature is also used to constrain the dust opticall depth along 
the line of sight. IRAS 17441-2411 shows no feature in IRAS LRS, we adopted 
amorphous carbon dust grains in the fitting. 

3. Results 

Table 1 lists all the fitting parameters we used. Figure 1 shows model fits to 
the SED for these three objects. Figure 2 shows the comparision between the 
observed model images. The model images not only reproduce the approximate 
shapes of the PPN, but also the absolute flux levels (as evidenced by the sizes 
of the outermost contour). In addition, searchlight beams in IRAS 17150-3224 
and dark lane in IRAS 17441-2441 are successfully reproduced as well. 



Table 1. Fitting Parameters 
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Figure 2. Comparision between the observed images and the model 
ones, (a) IRAS 17150-3224 (b) IRAS 18095+2704 (c) IRAS 17441- 
2411 All the grey-scale images are in log. scale, and the contours 
levels are the same in both observed and model images for each objects 
respectively. 



